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It was recently shown by computer simulations that a large class of liquids exhibits strong corre- 
lations in their thermal fluctuations of virial and potential energy [Pedersen et at, Phys. Rev. Lett. 
100, 015701 (2008)]. Among organic liquids the class of strongly correlating liquids includes van 
der Waals liquids, but excludes ionic and hydrogen-bonding liquids. The present note focuses on 
the density scaling of strongly correlating liquids, i.e., the fact their relaxation time r at different 
densities p and temperatures T collapses to a master curve according to the expression r oc F[p^/r] 
[Schr0der et al., arXiv:0803.2199 . We here show how to calculate the exponent 7 from bulk mod- 
ulus and specific heat data, either measured as functions of frequency in the metastable liquid or 
extrapolated from the glass and liquid phases to a common temperature (close to the glass transition 
temperature). Thus an exponent defined from the response to highly nonlinear parameter changes 
may be determined from linear response measurements. 



Introduction 



A liquid is termed "strongly correlating" if its viral {W) and potential energy ([/) equilibrium fluctuations correlate 
better than 90% 0, |3, [3 

at constant volume and temperature, 

AM^(t) = 7A[/(t). (1) 
where AW{t) = W{t) - {W), AU{t) = U{t) - (U). The "slope" is defined as 



and the correlation coefficient is defined as 



{AWAU) 

We have previously shown that the fluctuations of such liquids are well described by those generated by soft-sphere 
potentials (inverse power law potentials) [3, ^ . Strongly correlating liquids are approximate single-parameter liquids 
Ul) B El • Moreover the density scaling exponent is one third of the exponent of the approximate inverse power law 
potentials [sl. IT^. Recall that density (thermodynamic) scaling applies whenever the relaxation time t at different 
densities p and temperatures T collapse to a master curve according to the expression r cx F[p'^ /T]. It is generally 
reported now that density scaling applies for van der Waals liquids, but e.g. not for hydrogen-bonding liquids. This is 
consistent with our finding that the class of strongly correlating liquids includes van der Waals and metallic liquids, but 
excludes covalent, ionic, or hydrogen-bonding liquids - the latter three classes of liquids have competing interactions 
that spoil the WU correlation iS,, ilfl] . 

Fluctuation expressions 

Consider a viscous liquid with slow structural relaxation, i.e., with a relaxation time r that is much larger than 
one picosecond. The fluctuation-dissipation (FD) theorem for the frequency-dependent specific heat per unit volume, 
cv(lj), is given 3 by 

pOO 

keT^Vcviuj) = {{AEf) -iu {AE{Q)AE{t))e-''^'dt . (4) 

^0 

Relaxation takes place over a limited range of frequencies - typically 3-5 decades. By subtracting the responses at 
high ("00") and low ("0") frequencies well outside the relaxation frequency range, it follows that if tg is a time much 
shorter than those of the relaxations, but much longer than one picosecond, then 
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keT^V (cy(0) - cv(oo)) = {AE{0)AE{t,)) . (5) 

Write the energy as potential plus kinetic energy, E = U + K. Because the kinetic energy in the NVT ensemble 
fluctuates fast compared to one has {AE{0)AE{ts)) ^ {AU{0)AU{ts)). Thus, 

kBT^V (cy(0) - cy(oo)) = ((A[/)2)slow (6) 

where ((AC/)2)slow = {AU{0)AU{ts)) is the slow part of the potential energy fluctuations, i.e., slow compared to the 
picosecond time scale. 

Since the low-frequency limit gives the ordinary (dc) liquid specific heat and the high-frequency limit gives the 
"glassy" specific heat corresponding to perturbations that probe a frozen structure, this result may be written 

{{AU)Y°^ = fcBTV(4^^^'^ - 40lid) . (7) 
Similarly one finds for the virial fluctuations (where Kt is the isothermal bulk modulus) 

((AW^)2}slow ^ _ fcsTF(ifJ'l^^^ - /^folid) (8) 
and for the virial / potential energy correlation (where /3 is the pressure coefficient) 

Calculating R and 7 from data 

The WU correlation coefficient (Eq. [3]) on the tg timescale can now be expressed in terms of experimental linear 
response quantities as follows: 

, liquid _ ^solid^ 

R = ^^^^ ^ . (10) 

^_^^liquid _ ^solid)(,liquid _ ^soMyT 

Similarly the slope (Eq. ^ can be calculated as 
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_(^liquid_^S0lid) 

\ r(,uqnid_^sdid) • 



Calculating 7 for the commercial silicone oil DC704 

The slope 7 may be estimated from the high- and low- frequency limits of Cp{uj) and Ks{w) (and a rough estimate 
of ap) by proceeding as follows. Recall the identities 

cy = Cp - TalKr (12) 

and 



Combining these we get 



Kt = Ks — 

Cp 



(13) 
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TABLE I: 7 of DC704 data of from the high and low frequency Hmit of Cp{uj) and Ks{ui) af lid is estimated as a typical 

value of solids. 

~[K] 214 
4'^'^''^ [10^ J/(K m^)] IM 

4°li'l [10^ J/(K m^)] 1.05 
^hquid j^q9 3 g 

^|olid [^q9 5^ 

^liquid j^Q_3 j^-ij Q 5 

ll'l [IQ-^ K~'] 0.1" 

r(a|,^1^^d)24^1^^d/c^"''i 0.13 
r(aSolid)27^|olid/^solid q 

^hquid j^q6 j^3)j ^ 24 

^Splid j^q6 j/^j. jj^3-)] ^ q4 

^liquid j^q9 p^j 3^ 

j^solid [^q9 p^] 5.0 

7 (Eq. HDl 6 



(14) 



and 



/ TalKs \ 

KT = Ksil + ^^\ . (15) 

The value 7 = 6 is calculated from these expressions using our unpublished linear response data (table |T| . Unfortu- 
nately no experimental density scaling 7's are available for this liquid to compare to. 



OTP/OPP mixture and pure OTP 

The high- and low-frequency limits of the dynamic response can be estimated by extrapolation of static response 
functions of the glass and liquid phases to a temperature close to the glass transition temperature Tg. Table [Til lists 
extrapolated values of kt, Cp and ap for a mixture of o-terphenyl (OTP) and o-phenylphenol (OPP), Table Hill lists 
values for pure OTP. 

From the natural response functions of the constant pT ensemble it is straightforward to calculate the natural 
response functions of the constant VT ensemble: 

Kt = 1/kt, (16) 



Cv 



TalKr, 



(17) 



and 



I3v = apKr- 

Using these equations we arrive at the numbers in tables HIl and IIIII 



(18) 
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TABLE IL 7 and R of OTP-OPP calculated from data on Figure 4 of Ref. [H 

Tg [K] 233.7 

T/gjlQ-^niVmol] 203.9 

Cp'l''''^ [J/(K mol)] 364 

Cfoli'i [J/(K mol)] 236 
^Uquid p^_ij Q 

^solid [^0-9 p^-i] 

^liquid j^Q_3 j^-ij Q 74 

^solid K-i] 0.17 

^liquid jj^^j^ ^^j^j 284 

Cfpli'i [J/(K mol)] 229 
^Uquid j^q9 2.9 

^Solid j^q9 5 2 

^liquid j^q6 p^/j^j 2.1 
^solid [^q6 p^^^i q g 



R (Eq. [ini) 


0.8 


7 (Eq. [nj 


6.0 


^scale (j^gf ^) 


6.2 



TABLE III: 7 and 7? for OTP. Cp values are from Ref. is calculated from ap and Pv values from Ref. 

Tg [K] 244.5 
Vg [10~^mVmol] 206.1 



cJ^^[J/(K mol)] 336 

C|oli<i [J/(K mol)] 228 

^liquid j^Q_3 Q 7^ 

^solid [^0-3 K-i] 0.32 

^liquid j^q6 Pj^/j^j ^ 5g 

^^Olid |^q6 p^^^j 

^liqmd jj^^j^ ^^j^j 280 

CfPli'i [J/(K mol)] 209 

^liquid j^q9 p^] 2.2 

^solid [^q9 p^] 3 g 

7? (Eq. [To]) 0.3 

7 (Eq. [nil 4.1 

^scale ^^^^ 4 Q~ 



Summary 

We have shown that it is possible to calculate the density scaling exponent 7 from linear response measurements 
of specific heat and bulk modulus data. There are two ways to do this: Either by measuring broad-range frequency- 
dependent linear responses in the equilibrium metastable liquid phase or by extrapolations as done when evaluating 
the Prigogine-Defay ratio ^1^. Using the first method for DC704 we find 7 = 6. To the best of our knowledge there are 
yet no density scaling data for this liquid. Using the second method for the OTP-OPP mixture wc find 7 = 6.0 which 
compares favorably to the density scaling 7scale _ g_2 [l^; similarly we find 7 4.1 for pure OTP that compares 
favorably to the density scaling 7scale _ 4 q T^his good agreement may well be fortuitous given the uncertainties 
associated with our 7 estimates. Nevertheless these preliminary findings suggest that for strongly correlating liquids 
( "single-parameter liquids" ) the density scaling exponent - which refers to highly nonlinear parameter changes - may 
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be determined from linear response measurements. This is consistent with a general hypothesis of ours that strongly 
correlating liquids have simpler physics than liquids in general. 
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